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The possibility for detecting energy entanglement in normal metal-superconductor junctions is 
examined. For the first time we proved that two electrons in a NS structure originating from the 
same Cooper pair are entangled in the energy subspace. This work follows previous works where spin 
entanglement was studied in similar circuits. The device consists of a superconducting beam split- 
ter connected to two electronic Mach-Zehnder interferometers. In each arms of the interferometers, 
energies are filtered with coherent quantum dots. In contrast to previous studies of zero-frequency 
cross-correlations of electrical currents for this system, attention is drawn to finite-time measure- 
ments. This allows to observe two-particle interference for particles with different energies above 
and below Fermi level. Entanglement is first characterized via the concurrence for the two-particle 
spatial density matrix. Next, we formulate the Bell inequality, which is written in terms of finite- 
time noise correlators, and thereby we find a specific set of parameters for which entanglement can 
be detected. 
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I. INTRODUCTION 

Entanglement is the building block of quantum in- 
formation processing schemes In recent years, it 
has generated considerable excitement in the mesoscopic 
physics and nanophysics community. Over the years, 
electronic analogs of optical setups have been conceived 
in order to probe entanglement between electrons. In 
the context of electron transport, entanglement was first 
described theoretically 0, 0] in forked normal metal- 
superconductor devices revealing positive noise cross- 
correlations Q . Then entanglement in such devices was 
quantified via the Bell inequality test 

In all these schemes a superconductor injects Cooper 
pairs in two metallic arms, or in quantum dots connected 
to leads. When two electrons are emitted from the super- 
conductor, they either get into the same lead or in oppo- 
site leads. The resulting zero frequency noise correlations 
between the two normal arms will have a tendency to be 
negative in the former process due to the partition noise 
of Cooper pairs. The latter process, also called crossed- 
Andreev process, is the source of entanglement, and leads 
to positive noise correlations. Indeed, the Cooper pair is 
transmitted as a singlet state, and it is entangled both in 
spin and in energy. 

The detection of spin entanglement has been addressed 
on several occasions for ideal devices 0, Q , and with the 
assumption that parasitic effects are present @- In op- 
tics, a standard test for studying the degree of entangle- 
ment of two photons is to perform a Bcll/Clauscr Home 
inequality test, in a coincidence measurement. This test 




FIG. 1: a) Normal metal-superconducting fork with energy 
filters selecting electron (+E) on the left side and selecting 
holes (-E) on the right side, b) Same fork, but this time with 
spin filters, leaving the possibility of energy entanglement as 
both electrons and holes (with a definite spin for each side) 
can propagate in both branches. 



allows to probe the non-local nature of quantum mechan- 
ics: are the two particles in an entangled state? In con- 
trast, in condensed matter physics, experiments are often 
performed in a stationary state, with constant currents 
flowing at the output. Nevertheless, in Ref. Q it was 
shown that a Bell inequality test [!j can be expressed 
in terms of zero frequency noise cross-correlations at the 
output. Subsequently, similar Bell inequality tests were 
proposed for orbital entanglement in a superconducting 
setup m normal metal devices flTL f]~3 | . and in Quan- 
tum Hall systems HE EH- 

Ref. proposed that some constraints should be im- 
posed on a plain superconducting-normal metal fork Q in 
order to probe entanglement. Energy filters where added 
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on each side [Fig. ma)], selecting positive/negative ener- 
gies depending on the side. An incoming hole thus could 
only be reflected as an electron in the opposite branch, 
yet the spin degree of freedom was untouched. Energy 
entanglement was also proposed Q using spin filters in- 
stead of energy filters [Fig. IHb)]. Except for the charac- 
terization of positive noise correlations which constitutes 
a symptom of entanglement, so far no quantitative test 
of this entanglement has been reached. 

Recent advances in quantum optic experiments on 
momentum-phase entanglement |l5j led to experimental 
realization of such techniques as two-particle teleporta- 
tion |l6j ] , t he p urification protocol for polarization entan- 
gled state |l7j and the full Bell-state analysis 0] ■ All of 
the above provide additional motivation for us to reex- 
amine the question of energy entanglement. 

The purpose of this work is to explain how en- 
ergy entanglement can be detected in normal metal- 
superconducting forks. Although the setup which we 
propose is built from elementary building blocks of meso- 
scopic physics, it is rather elaborate. We thus emphasize 
that our main goal is to show whether energy entangle- 
ment detection is plausible from a theoretical standpoint. 

The paper is organized as follow: Sec. 2 gives a basic 
description of the setup, which is composed of a supercon- 
ducting fork connected to two Mach-Zchndcr (MZ) inter- 
ferometers. Sec. 3 and 5 describe the calculation of the 
two-particle density matrix and electrical current cross- 
correlations in the context of scattering theory for this 
normal metal-superconductor setup. Sec. 4 describes 
the entanglement characterization via the concurrence of 
the density matrix. Sec. 6 develops the Bell inequality 
analysis in terms of electrical currents. 



II. DESCRIPTION OF THE SETUP 

We now present a setup (Fig. |5Jl, which is able to test 
energy entanglement via the violation of Bell inequali- 
ties. We use a normal-metal fork connected to a super- 
conductor. The leads connected to the arms of this fork 
are assumed to have quasiparticle mean free path larger 
than the size of the setup: in this ideal version of the 
setup, backscattering is absent. A possibility would be 
to use electrostatically defined quantum wires on semi- 
conductor heterostructures, assuming that the interface 
between the two dimensional electron gas and the super- 
conductor can be controlled. The two leads connected 
to the injection fork (center of Fig. |2J) are each attached 
to two other forks (left and right of Fig. |2J) , which filter 
the particles according to their positive/negative energy, 
measured with respect to the chemical potential of the 
superconductor. This energy filtering can in principle be 
realized via a double barrier structure forming an elec- 
tronic Fabry-Perot interferometer. An example of such a 
filter was demonstrated in a "which pat h detection" ex- 
periment at the Weizmann Institute |l9j ^ where one arm 
of the devices contained a coherent quantum dot. 




FIG. 2: Setup for measuring energy entanglement from a nor- 
mal metal-superconducting fork. The superconductor emits 
Cooper pairs in the ballistic region. The electrons from a pair 
are then spit on the first fork. Then the particles are filtered 
according to their energy. The amplitudes for particles with 
energies above and below Fermi level are combined at the last 
beam splitter, in close analogy with the optical setup |T^I . 



Next, the two forks on the right and on the left are 
reconnected, forming MZ interferometers. The paths on 
each arm of such interferometers have equal length, and 
the two beam splitters are assumed to be symmetrical. 
Later on in this paper, we shall consider the measure- 
ments of current cross-correlations performed between 
different contacts A\, A2, -B3 and B4 of these two in- 
terferometers. Note that, in this suggested geometry of 
the setup, there is no absolute need for spin-filters, con- 
trarily to what was suggested in Ref. and illustrated in 
Fig- Gib). As we shall see below, if the interface between 
the injecting fork and the superconductor is opaque, the 
two constituent electrons of a Cooper pair are split be- 
tween the right and left side and the cross correlation 
signal is positive, as in a photon coincidence measure- 
ment. 

The idea of this experiment is that two electrons, which 
originate from a Cooper pair in the superconductor, once 
on the normal metal side, have energies which are sym- 
metric with respect to the chemical potential of the su- 
perconductor. In principle they can be filtered in a way 
that an electron with energy above this Fermi level always 
gets into the upper arms of the interferometers (energy 
+eo, leads ua and ub in Fig. |5J) or it gets into the lower 
arms of interferometers otherwise (energy — eg, leads (La 
and ds in Fig. |2J). That is, two electrons being split at 
the injection superconducting fork (center of Fig. |2J) al- 
ways end up into the right upper and the left lower arms 
of the interferometers or into the right lower and the left 
upper ones. 

This situation is very similar to the phase-momentum 
entanglement in quantum optics. Indeed, an experimen- 
tal non-locality test for photons was achieved by Rarity 
and Tapster using two MZ interferometers |15| . Two pho- 
tons with different wave lengths were generated during 
a down-conversion process, and subsequent recombina- 
tion of the beams allowed to perform a correlation mea- 
surement. We find that this quantum optics experiment 
can be translated for and applied to electronic devices. 
As we shall sec, in our condensed matter analog, which 
uses charged particles, the Bell inequality violation can 
be tuned by varying the magnetic fluxes within intcrfer- 
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ometers (electrons gain extra phases e ±ltPa / 2 and e ±l(pb ^ 2 
in the lower (upper) arms of the interferometers). 

In the quantum optics experiment [lfij . the momen- 
tum direction of the down-converted photons allowed to 
separate the photons and to let them recombine at the 
beam splitters, provided that a single beam splitter can 
be reached only by photons with equal frequencies. This 
property is essential for the two-particle interference. In 
our proposal the particles, which are reaching the last 
beam splitters (recombination of the beams) along dif- 
ferent paths, have different energies, so it is in principle 
possible to achieve "which path" information and there- 
fore to destroy the interference. We propose to use mea- 
surements on short enough time scales so that it is not 
possible to resolve the energy of the particle (energy/time 
Heisenbcrg principle) . We claim that the two-particle in- 
terference and the non-local properties of quantum me- 
chanics can be detected in this ideal setup. 

III. CALCULATION OF THE DENSITY 
MATRIX 



The states Uij/3 or Vijp correspond to wave functions of 
an electron or a hole scattered in lead i, due to a quasi- 
particlc of type (3 = e,h incoming from lead j. Opera- 
tors Cjf3 a (E) and c^^E) satisfy standard Fermi statistic 

relations. v{{E) — hk{(E) and v J h {E) = hk J h (E) are ve- 
locities of electrons and holes in lead j. 

The evolution of particle states Uijp and Vijp is de- 
scribed by Bololiubov-de Gennes equations, which can 
be written as 



f ti 2 d 2 \ 

+ A(x)vijp(x) = Ev,ijp(x), (3a) 
/ h 2 d 2 \ 

+ A*(x)u,ijp(x) =Evijp{x). (3b) 



In order to analyze the properties of the injected elec- 
trons, we calculate the two-particle density matrix for 
quasiparticles caught right after energy filters (in leads 
uai d-A and ub, ds)'- 

g(a,x) = Tr [p^x^l^x^ a2 {x 3 )i> a ^)\ , (1) 

here we have written the whole matrix with non-diagonal 
elements, and the indexes a describe the spins part of 
density matrix, and orbital indexes i = ua,(Ia and 
j = UB,ds describe the lead there the particle is caught. 
In practice, only diagonal elements with Xi = x^ and 
Xj = Xjt can be measured, but our task is to analyze the 
amount of entanglement which is implicit in such a ma- 
trix, therefore non-diagonal elements play a crucial role. 

In our calculations we follow the scattering theory ap- 
proach developed in [2(], El, an d later applied to NS- 
systems with Andreev scattering j2^. For background 
purposes, the reader is referred to the reviews pH l24j . 
Here we follow closely the notations of Ref. |2^. For sim- 
plicity, we consider ballistic quasi one-dimensional quan- 
tum wires without backscattering. The transport in this 
system is governed by the properties of Andreev reflection 
and normal reflection of electrons and holes. So it is use- 
ful to perform Bogoliubov transformation (see Ref. 126^ 
for the annihilation operator of a particle at a position x 
with up/down spin a = ±1: 




The pair potential A(x) should be calculated selfcon- 
sistently, but, for simplicity, in our calculation it cor- 
responds to the superconducting gap A in the bulk of 
the superconductor (x > 0), and it is zero in the normal 
leads (x < 0). 

In a normal ideal lead A(x) = and V(x) = 0, 
Bogoliubov-de Gennes equations l|3a[l - l|3bjl reduce to 
Schrodinger equation for free electrons and holes. So- 

for holes 



Cjj3a{E, ti)- 



lutions of form e lk ? x for electrons and e 
are chosen, where 



\/2m(/ic 



E)/U and k% 



(2) 



y^2m(ps ~ E) /h are wave vectors for electrons and holes. 

Electrons and holes that originate from a particle of 
type (3 in lead j and scattered into lead i are described 
by: 



Uij0 ( Xi ) = 6, tJ 5 ef3 e + ^^ + StfJ -^e-^*, (4a) 
v ijP ( Xi ) = Mwe-<*?*< + StfJ (4b) 



Here we have the opposite sign for momenta of electrons 
and holes, which is due to the symmetry of Bogoliubov- 
de Gennes equations. For simplicity, we introduced a no- 
tation S^J 3 for scattering-matrix element expressing the 
amplitude of an outgoing particle of type a in lead i due 
to an incident particle of type (3 in lead j. In general, 
these amplitudes depend on the quasiparticle energies. 
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We linearize the wave vectors in energy in the exponents, and then take them equal to Uf in pre-exponential factors. 



We substitute all the equations into Eq. QJ and obtain: 

(4<t; ( E Wti')^ a > 3 (E2,tj<)Cmo- 2 (E3, tj)c nai (£4, U)^ U*, k (Xi>)Uj, l (Xj> )u jm (Xj)Ui n (Xi) (6) 

+ \clv> 1 ( E uU')ci-<T' 2 (E2,tf)ct n _„ 2 (E 3 ,t^ (7) 

+ ( 6 L; ( E i^i')ci-a' 2 (E 2 ,t jl )c m(T2 (E 3 , t 3 )cl_ ai (Ei, t*)) ^i^'k^^i'iix^u^x^vl^Xi) (8) 

+ (c/c- CT ; {E u t v )c\^(E2, V) 6 L-a 2 {E? n t : j)c n!yi (Ei, i;)) a[a 2 v vk (x v )u* jn (xj>)v* jm (xj)ui n (xi) (9) 

+ (^Ck-a[{Ei,t i /)cl a , 2 (E 2 ,t jl )c mcr2 (E 3 ,t j )cl l _ (7i (E 4 ,t l )^ a'^iVi'kixi^up^Xj^Uj^x^v^Xi) (10) 



+ ^4-^(^1, iiOci-^^^jOC-^^^jOcl-o-i^,^)) o-i^t7 2 o-iUi'fc(a: i ')?; 3 v i (a; J -/)?;* m (xjX„(a;i) , (H) 

here the summation over indexes k, I, m and n implies a summation over lead indexes and types of particles - e or h. 

We consider the case where all normal leads are biased with the same voltage V. With these assumptions, the 
density matrix takes the form (see Appendix): 

g((J,x) — d{j 1 (7 / 1 d(7 2 (T 2 Gi'iGjj' $0ia' 2 $a 2 (T\Gi' "i~ ^CTiaz^-tr'^CT^Ei' j' E^j , (12) 

where we have used a notation for a matrix: I ai a 2 = ify and fj, is a Pauli matrix. We compute this density matrix at 
locations taken in the leads of the interferometers: ua, cIa, ub and ds (See Fig- - We thus define indices i and i' 
which can take the values ua or d,A and indexes j and j' equal to ub or ds- We also use a notation for the two-point 
correlation functions: 

1 f + °° 

C, ,[.<■•■■'■. - ! j - t i) = l—y\ / ^/fc(^)«fe(^)«^(^))e" if!( * 3 "* l) , (13a) 

nVF ~ J-oo 

1 /•+<» 

F tJ { Xl ,x 3 ,t 3 - U) = t—Y] / ^/ fe (£)« fe (^)^ fe (^)) e - tB ^- t '). (13b) 

^ J_ OQ 



Up to now, we have not specified the trans- 
mission properties of the combined normal metal 
fork/superconductor plus energy filtering forks. For the 
parametrization of the beam splitter connected to the su- 
perconductor we take a scattering matrix with real am- 
plitudes described in |27j : 

f-(a + b) y/i yfe\ 
U SLo r = a b ], (14) 

V v 7 ^ b a j 

with the free parameter < e < \ describing the trans- 
mission probability between the top lead and the right 
and left arms. In general, there are two possible choices of 
parametrization Q : a = n/1 ~ 2 2<; ~ 1 , b = v/1 ~ 2 2e+1 , or alter- 
natively a and b can be exchanged. Here we use the first 
parametrization: it corresponds to a good transmission 
for electrons (holes) entering from the left lead of the fork 
and exiting from the right one, and vice versa. We claim 
that the excess parts of all the correlators are the same if 
one choses the opposite convention. This unitary matrix 
(|14|l is assumed to be the same for electrons and holes, 



r 



which implies a weak dependence of the scattering matrix 
on energy. Due to multiple Andreev reflection on NS- 
boundary and a phase shift in the superconducting lead 
(S), the transmission amplitudes change in an analogous 
manner to a Fabry-Perot interferometer. Assuming ideal 
Andreev reflection on the boundary of the superconduc- 
tor: rf s = r h s e s = e - iarccos # , and e _i7? = rf s e~ iEr ^ 



one can get 



e(n + b)p~' 2irl 

A = rgg=rg% = a- K , } * 9 . , (15a) 
LL HH 1 - (a + 6) 2 e~ 2 "' 

d _ fata _ f aa — l e(a + b)e 2 "' 

B-t LR -t RL -b 1 _ {a + b)2e _ 2iv , (15b) 

C = rf L = r RR = 1 _ {a + b)2e ^ n > ( 15c ) 



D = tf R = tf r = , , N9 9 . , (15d) 

LR RL 1 _ ( a + b) 2 e- 2tr > 



here A and B arc reflection and transmission amplitudes 
for particles of the same type (aa) with incident and fi- 
nal lead to the left (L) and to the right (R) from the 
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fork, and C and D are reflection and transmission am- 
plitudes for the case of changed types of particles (a/3). 
Although notations A, B, C, D are not used here, they 
arc introduced for later convenience, in Eqs. H35a|) and 
(l35bll . 

The energy filtering forks are chosen "ideal" : electrons 
with energies above the superconducting chemical poten- 
tial pass without reflection through the upper arm, while 
electrons which energy is below this chemical potential 
go in the lower arm of the fork (left side of Fig- El - The 
following transmission amplitudes are thus chosen for the 
left energy filtering fork: 



Ax,- 



Xi and Atij = t j 



tf, and 5ij (Sf-) is equal 



rf L {±E) = 0, 
rZ(+E) = +1, 



,(- 



t e u l(+E) 



1- 

0, 



E) = - 
r%(-E) = 0, 
C L (-E)=0, 

t e u L ( 



1. 



-E) 



(16a) 
(16b) 
(16c) 
(16d) 
(16e) 



here the amplitudes for hole type quasiparticlcs follow 
from Eq. (|A16|I and the same is correct for the right fork 
exchanging L on R. This is a very simple parametriza- 
tion, and the only constraint is the unitary conditions. 
In a real experiment there may be a more complicated 
dependence on energy leading to imperfections, but this 
question goes beyond the scope of the paper. 

Now the global scattering amplitudes are calculated in 
a simple way, because of our assumption of no backscat- 
tcring in the leads. For example: 



1L LL Lt 



These amplitudes allow to calculate the pair correla- 
tion functions from the Eqs. (|13a[) and l|13b[) . and the 
results may be divided into an equilibrium and an excess 
part: 

Gij (xi , Xj , At^ ) = G e t j (xfj , Atij ) + G^j (Axij , At io ) , 

(18) 

the former depends on the sum of coordinates and de- 
scribes the Fermi correlations within normal leads, the 
latter depends on the difference of coordinates and de- 
scribes the particles injected from the superconductor. 
The most interesting contributions for the discussion be- 
low are the excess contributions which dominate at low 
temperatures (T <C eV): 



ij 



zpex 



'2(1-6)2 



n(Axij + AUj) 



6Vl-2e +ikFX f ; JH^tA^+Att 
% 2 (l-e) 2 

eV 



n(Axij + At^) 



(19a) 



(19b) 



where i = ua^a and i = ub^cIb', e is the scatter- 
ing parameter of Ea. (|14f> of the superconducting fork; 



to unity (zero) only if i and j are simultaneously u or d 



IV. CONCURRENCE FOR THE DENSITY 
MATRIX 

We now turn to the calculation of the concurrence, as 
defined in the work of Wootters [2^| , but adapted to our 
transport geometry. With the additional assumption of 
simultaneous measurements (all the times ti in Eq. (|19bl) 
are the same) and equal coordinates on the same sides 
of the setup (x UA = Xd A and x UB — x c i B , up to a the 
precision Sx <C vpTy) the whole density matrix simplifies 
to: 



g(a,x)= g% [S^S^Ji'iSjf 

-<Wi<W; s i'j$i3' e±tip g 2 ( A x) 
+I^ 2 I^ 2 5 t , 7 6 i] e ± ^f 2 (A x ) 



where 
3o 



eV 



2(1 -e) 2 2-Kvph 



, ff(Aar) 



sm 



Ax 
2v f t v 



Ax 



/(Ax) = y-LJlg(Ax), 



(17) where Ty = h/eV. For i 



J 



u and i' 



J 



(20) 



(21) 



(22) 



d, the 



parameter <p 



Ax which appear in Eq. ||2B is taken 



with a plus sign, for i = j' = d and i' = j = u - with a 
minus sign, and in all other cases it is zero. 

The two-particle density matrix describing the spin en- 
tangled state of two quasi-particles inside a superconduc- 
tor has been already studied in |29j | for a bulk supercon- 
ductor. In contrast, our analysis allows to study energy 
entanglement in a condensed matter setting with nor- 
mal metal leads, which are essential for the diagnosis of 
entanglement. Since we are not interested in the spin en- 
tanglement here, we trace Eq. (|20|) over spin degrees of 
freedom, and we obtain a density matrix which depends 
only on orbital degrees of freedom: 



gi,j-,i',j< = Tr a g(a,x) 



2glNp AB 



(23) 



Pab = -77 











2 + / 2 

(f 2 -g 2 )e^ 






{f 2_ g 2 )e+iv 

2 + / 2 








2-.9 2 
(24) 

where the normalization factor N = 8 + 2/ 2 — 2g 2 is cho- 
sen so that the trace for this matrix is unity. The param- 
eter g 2 < 1. It turns out that the parameter ip = -^-Ax 
does not appear in the final result for the entanglement 
measure. 
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One of the possible expansions (on the basis of pure 
states) for such a mixed state may be written in the space 
of the pseudo-spin | f) = \u) and | J.) = \d): 



1 

PAB = jj 



(/ 2 -.9 2 )l* (+) )(* (+) l + (2-ff 2 )i 



+2g 2 (| UIb){UIb I + I UTbXUTb I)] , (25) 

where I - the unity matrix, and the only entangled term 
in the sum corresponds to the state: 



I* 



(+)\ 



1 

71 



(I UW + e^l UTs)) 



(26) 



Note that the density matrix of Eq. (|TJ describes only 
a part of the whole wave function for all electrons, and 
this is why it corresponds to a mixed state. But for the 
case of an opaque injection fork (e — ► 0, f/g — » oo), the 
state of the orbital degrees of freedom is represented by 
the wave function (|26() . which corresponds to one of the 
maximally entangled states. Thus, according to the so- 
called " monogamy of entanglement" j^Ol ; these degrees of 
freedom cannot be entangled to any other ones, and one 
can say that these particular two electrons are maximally 
entangled. 

We now calculate the concurrence for pab - an entan- 
glement measure proposed by Wootters et al. pgj : 



C = max{0, v Ai 



A 4 }, (27) 



where Ai > A 2 > A 3 > A 4 > are the eigenvalues of 
the matrix pabPabi an d the " spin- flipped" matrix is the 
folloing: 



- P AB=t { y A) ®^ ) P\B* { y A) 



T (B) 



So for our setup the concurrence takes the form: 

f-2 



C = max < 0, 



r- 



or using the ratio between / and g: 

(l+ 7 ) 3 2 -2 7 



C = max < 0, 



g 2 + 4 7 



(28) 



(29) 



(30) 



where 7 = 1 _ 2 e e _ e2 is the transparency parameter of the 
superconducting fork. 

The concurrence, which depends on the coordinate dif- 
ference , is shown in Fig. [3] choosing 7 = 0.001. It 
takes values between and 1, the former corresponds 
to absence of entanglement, and the latter to maxi- 
mal entanglement. Actually, for 7 going to zero and 
Ax <C vftv the concurrence is maximal, which means 
that all the excess cross-correlation noise is due to the 
pairs of electrons in |^( + )) maximally entangled orbital 
state. 

A natural question to ask is whether entanglement can 
occur between orbital and spin degrees of freedom. Start- 
ing from the general density matrix of Eq. I|20|l . it is 




FIG. 3: (Color online) Plot of the maximal Bell inequality 
violation B max (middle green plot) as a function of the de- 
viation from the coincidence of the current measurements: 
a = XB ~ XA for the case of small measurement times: r <C 
h/eV. B max is sandwiched between the plot for concurrence 
C (lower red plot) for orbital part of the two-particle density 



matrix: pab, and the plot for the function: 



(upper 



blue plot). The doted line at l/\/2 marks the critical value 
above which the Bell inequality is violated. The transmission 
probability is taken to be: 7 ~ e 2 = 0.001. 



possible to perform a partial trace with respect to a sin- 
gle spin and a single pseudo spin, in order to test this 
property. We do not show this computation, but instead 
mention only the result: the concurrence is strictly zero. 
In any case, if no trace of the density matrix 1)20(1 is 
performed, we conjecture that the joined state of spin 
and orbital degrees of freedom can be a hyperentangled 
state |3lj . In quantum optics these states can be used 
for instance to measure the concurrence directly |32j| . 



CALCULATION OF CURRENT 
CROSS-CORRELATIONS 



The theoretical calculation of the concurrence cannot 
easily be tested experimentally: it would require quan- 
tum state tomography, which represents a considerable 
challenge in experimental mesoscopic physics. A theo- 
retical proposal for tomography has been been presented 
recently [33j . which relies on the measurement of zero 
frequency noise cross-correlations. Because of the con- 
siderable experimental challenge of tomography and be- 
cause we wish to compare different diagnosis for quantum 
mechanical non-locality, here we present a Bell violation 
test based on the analysis of short time current cross- 
correlations. Here we are interested in the same quanti- 
ties for the specific setup of Fig. [21 evaluated at points in 
the leads: A\, A2, B3 and B4. 

We continue using scattering theory for NS-systems. 
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Our starting point is the current operator in lead i: 



dx 



dx 



(31) 

this operator is expressed in terms of electron and hole 
wave functions: 



Ii(x) 



ieh 1 



2mvp 2 



i p + oo p + oo 



[( 



~ (J \ v in Cf x u irn ~ Cf x v in U im) C m(T C n~a 



(32) 



where we have replaced the sums over j and (3 by a single 
sum other index m. The terms with index m (n) corre- 
sponds to energy E\ (E2). The chemical potential fis of 
the superconductor is large compared to any other energy 
scale in the system, thus the assumption k e = kh = kp 
is made here, which explains the presence of the Fermi 
velocity Vp in Eq. 

The calculation of current correlations is rather stan- 
dard, but it is reproduced in the Appendix so that the 
notations of the present work can be understood. Com- 
pared to the calculation of the density matrix of Sec. 3, 
the characterization of the scattering properties of the 
setup of Fig. |2 now requires to include the MZ interfer- 
ometers. The magnetic fluxes which are applied on the 
left and on the right ($_b) interferometers define 

the Bell angles: 



2n- 



Vb = 2tt^-, (33) 



In analogy with spin entanglement of electrons 0, , 
where electrons are detected via spin filters with arbi- 
trary polarization, the phase angles if A and Vb allow to 
rotate the pseudo-spins with respect to the Z axis, which 
is defined in a basis | t) = \u) and | J.) = \d). Indeed, 
electron wave functions gain a factor e ±lVA - B ^ 2 in the 
upper/lower arms of the interferometers, which is equiv- 
alent to a rotation by the angle <pa,b in the pseudo-spin 
basis. Next, the use of a symmetric beam splitter, which 
is described by a unitary matrix, 



U= f ± 

\ V2 V2 



(34) 



performs a rotation of such pseudo-spin around the Y- 
axis by an angle 7r/2. In the end, the latter measure- 
ments at the points A\,Ai (left side of Fig. EJ) and -E>3,£?4 
(right side of Fig. 0) correspond to a Stern-Gerlach type 
of measurement along Z axis for the initial pseudo-spins. 
These transformations on both the right and the left or- 
bital pseudo-spins allow to perform standard Bell type 
correlation measurements. 



For constructing the scattering amplitudes, we use A, 
B, C, D from Eqs. 115a|) - (|15d(l and the elements of the 
beam splitter scattering matrix Eq. I)34|l. The follow- 
ing convention is used for the angle signs: an electron- 
like quasiparticle on the left (right) side gains a phase 
e +iVA/2 ^ e -up B /2^ wncn propagates in the clockwise 
direction through the upper arm of the interferometer. 
For anticlockwise movement the sign in the exponent is 
opposite. For a hole-like quasiparticle moving clockwise 
in the lower arm the sign is also opposite. 

Note that this convention applies to the case where the 
direction of the magnetic field is opposite in different in- 
terferometers. This situation may seem hard to realize 
experimentally, however it was chosen only for conve- 
nience: in the contrary situation (with both fields in the 
same direction) one simply needs to change the sign of 
an angle in 

For a given type of particle (electron or hole), the scat- 
tering matrices have dimensions of 4 x 4, which corre- 
sponds to the number of leads attached to the supercon- 
ductor: recall that here we are working in the Andrcev 
regime. In order to write the matrices, we use the ampli- 
tudes in Eqs. I|15a|) - (|15d|) . an angle Sip = tpA — <Pb and 
unitary matrices for scattering amplitudes of the beam 
splitters in Eq. (|3lJ) : 



/ 



S aa {+E) = 



for a 7^ [3: 



S a/3 (+E) = - 



2 ~ 2 

iA i_ 

B 2 +iS?/2 
2 c 

iB +iSy/2 
2 e 



iA i_ 

2 2 
__A _ 1 

iB^+iS<p 2 /2 
2 e 
_B +iS v /2 
2 e 



i8ip/2 
iS v /2 
1 



2 ~ 2 
iA i 
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c 
1 

iC_ 

D iS V /2 
2 c 

iD iSip/2 
2 C 



2 

_c 

2 

iP c iSip/2 

_ 2 P iSip/2 
2 c 



iSip/2 
-i5ip/2 



C 
2 



iC 
2 



i_B_ -iSip/2 
2 e 

B -iS<p/2 
' 2 c 
iA i 
2 2 
A _ 1 

(35a) 



,-i<5^/2 

e -iS v /2 

iC 
2 
C 

2 (35b) 



In these matrices, the terms ±^ and +| in the scattering 
amplitudes originate from the full reflection of a particle 
approaching an ideal energy filter, as specified by Eqs. 
(|16b(l - H16cfl . which is adjusted not to allow the transmis- 
sion of a particle with this particular energy. The ampli- 
tudes for energies below Fermi level (— E) are calculated 
using the property of Eq. (|A16|) . 

Note that here we did not include explicitly factors 
such as e ±k '/ hl , which correspond to phases accumulated 
by particles propagating through a lead of length I. These 
phases can be included as an overall multiplication fac- 
tor on the scattering matrix in Eqs. l|B5a(l - (|B5c|l . de- 
pending on the lead coordinate. Here the assumption 
that both arms of any given MZ interferometer have the 
same length is implicit. We also neglect the length S 
of the injecting lead of the superconducting fork, and we 
assume an ideal NS-boundary in the pure Andreev reflec- 
tion regime with amplitudes: 



1 ss 



ss 



From the scattering amplitudes of Eqs. (|35a|l - lj35b|) we 
obtain the noise cross-correlations, which may be sepa- 
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rated in an excess and an equilibrium parts: 

C{Xi , j Xj — C ex {x% , ti , Xj , tj ) H~ Ceg (-^i s ^2 5 -Ej 7 )) 

the excess part (first term) depends on voltage and orig- 
inates from the first term in Eqs. (|B5a() - (|B5c(l . the equi- 
librium part (second term) is voltage independent and it 
occurs from the second term (delta functions) in Eqs. 
(lB5aV(|B5cj) . 

It is then natural to introduce time scales which char- 
acterize the time spacing between electron wave pack- 
ets, as well as the time of flight of particle through 



the wires: ry = h/eV, t { . 



i\x^\ + \xi\)/v F and 



j\ — \xi\)/vF- In addition, r,j = tj — ti. So 
the cross-correlations at finite temperature T reads: 



( T ij > T ii ; T) 



2e 2 e 2 (l-2e-e 2 ) 



v VlJ1 13 h 2 8(1 -e) 4 

x(l + (-l)^cos(^-(r y -r-)/r y )) 

(Tij ~ Tff) (ttT) 2 



x 4 sin 



2t v sinh 2 (7rT(ry -T??)/h) 



,(36) 



2e 2 1 + VI - 2e - (2 + y/i - 2e)t 



ij \'W> 'ij > ^ I ^2 

(ttT) 2 



:i- £ ) 2 



sinh 2 (7rT(Ty + r,f/)/ft) sinh 2 (7rT(r lJ - t^)/K) 



For zero temperature (fcsT <C eV), the correlators re- 
duce to: 



2e 4 V 2 e 2 (l-2e-e 2 ) 



2\ sin 



V 2tv 



/i 2 



;i- £ ) 4 



2t v 3 



x(l + (-l)' i+ ^'cos (ty-fa. 



ex 



)Av)), (37) 



VI. BELL TEST FOR CURRENT 
CORRELATIONS 

In optics, the Bell inequality test is typically expressed 
in terms of correlators of number operators. This re- 
flects the fact that a coincidence measurement on two 
photons is performed. In Refs. and ||, these num- 
ber correlators were expressed in terms of current noise 
cross-correlators, which contain an irreducible contribu- 
tion and a reducible contribution. The reducible con- 
tribution is proportional to the product of the average 
current in each lead. The irreducible contribution domi- 
nates for short observation times. Subsequently, the Bell 
inequality can be written in terms of irreducible zero fre- 
quency noise cross-correlators between the different arms 
of the device. A slightly different point of view was pro- 
posed in and discussed in : As long as the spacing 
in time between electron wave packets e/(J) is large com- 
pared to the size of these wave packets H/eV, successive 
pairs do not mix so that all information about entangle- 
ment is included in the zero frequency noise correlators. 

In this letter, we suggest a Bell type inequality for finite 
time cross-correlations of electrical currents. It is similar 
in spirit as in work of Ref . 0, 0, where Bell inequalities 
are directly expressed in terms of current correlators, yet 
we shall see that our analysis is not necessarily restricted 
to short times. We thus follow the derivation of Ref. 
for the Bell inequality violation with particle numbers: 

\E(tp a ,tp b ) + E{ip a , <p' b ) + E( 9 ' a , 9b ) - E(tp' a , <f/ b )\ < 2, 

(39) 

where 

w( v ((iVi(T)-A 2 (r))(iV 3 (r)-iV4(r))) 
^'^ )= ((A. l( r) + A 2 (r)) ( A 3 (r) + A 4 (r)))' (4 ° } 

and with ip a , ip b - angles corresponding to different mag- 
netic fluxes in the left and the right interferometers, r is 
the duration of the measurement in time. Particle num- 
ber operators in different leads A1-B4 arc defined from 
current operators: Ni(t,r) = J* +T dt'Ii(t'). We perform 



CSW-C 0) = - 



2e 2 1 + \fl — 2e — (2 + \fl — 2e)e density matrix averaging and time averaging: 



(27T) 2 



[l-ef 



r eq\2 



(38) 



+T 



Note that the equilibrium noise describes Fermi cor- 
relations within the normal leads. At low temperatures 
it decays as an inverse square with large distances from 
the superconducting fork to the detectors. So assuming 
large distances, Xi, Xj 3> Wftv, only excess noise is rel- 
evant, and this corresponds to the excess noise due to 
Cooper pairs injected from the superconductor. The de- 
pendence on magnetic fluxes in the excess noise occurs 
via Sip = ip a — ^Pb- One should therefore attempt to 
construct a Bell type inequality with these measurable 
quantities. 



(N^N^t)) = lim — / dt(N t (t,r)Nj(t,r)) . 

T— >0O Zl J _ rp ^ 

(41) 

We define the irreducible current cross-correlator at 
given positions and times t\ and t 2 : 



Cij(t 2 -ti,Xi,Xj) = ^Ii(Xi,ti)Ij(Xj,t2) 

^(x,,t 1 ))(/,(x J ,t 2 )), (42) 

and we write the Bell term in the following form: 

£ dh £ dt 2 (C 13 - C 14 - C23 + C24) + A- 



E(lfi a ,lfi b ) 



/ T dti J T dt 2 (C 13 + C14 + C 23 + C* 24 ) + A+ 



(43) 
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where A± originate from the reducible part, which is typ- 
ically proportional to the average currents: 



A- 



((h)±(I 2 ))((I 3 )±(h)) 



(44) 



The coordinates Xi and Xj in Eq. 1|42|) correspond to po- 
sitions of the detectors counted from the superconducting 
fork, and r correspond to duration of measurements. 

For the case of symmetrical beam splitters used in in- 
terferometers and a symmetrical superconducting fork 
the average electrical currents are identical: (Jj) = 

^TT" (i-e) 2 ; so A_ = and 



A+ = At 



2 eV 



h {l-e) 2 ' 
In this case the Bell term reads: 



(45) 



E(ip a ,ipb) = 



in 2 ^-H + c 
(t 2 -t!+Q) 2 



cos (Sip — (t2 —ti+a)) 



(t 2 -t 1 +a) 2 



+ 7/32 



here a = XB x 

V F T V 

where x\ — x 2 



(46) 

For simplicity, we consider the case 
xa and x% = X4 = Xb, assuming a 
precision Sx <C vpTy - i.e. the detectors of each party 
A and B are equidistant from the superconducting fork. 
Moreover, (3 = t/t v , 7 = 1 _ 2 £ e _ e2 . Sip = ip a - ip b . So 
there are three dimensionless parameters: a - which de- 
scribes the fact that the measurements performed by A 
and B parties arc not in coincidence; (3 - the duration 
of these measurements, 7 - the parameter describing the 
transport properties of the superconducting fork (recall 
that it depends only on e, which is the probability of a 
particle to go from the injection lead into one of the two 
leads); and Sip is the difference between the Bell angles 
introduced in Eq. (|33() . 

In our analysis of the Bell inequality of Eq. (|39[) we 
consider two simple cases, for both of them the Bell term 
in Eq. l|4"U|) satisfies E(ip a ,ipb) ~ cos (ip a — ipt). 

i) The case of coincident measurements: a — 0. On 
Fig. 0] we show the resulting dependence of the maximal 
Bell inequality violation B max (which is the maximum 
of the left side of the Bell inequality over the Bell an- 
gles) over the measurement time r. From Fig. 0] it is 
clear, that for the case of superconducting fork with very 
bad transmission (e — > 0, 7 — > 0), the violation of Bell 
inequality is possible only for r < 2ry. We again stress 
that long measurement times destroy the interference be- 
tween trajectories with different energies, so that the Bell 
inequality cannot be violated. 

In our analysis, the maximal violation depends on the 
transmission properties of the superconducting fork: the 
closer 7 to zero the better the violation. The Bell in- 
equality in Eq. (|39[1 may be violated maximally (with 
(3 — > 0) for e — > 0, and it is not violated at all for 
e > tmax = (3 — \/2)/7 ss 0.22 [this corresponds to: 
7 > Imax = (V2 — l)/4]. This fact may by explained in 
the following way: the irreducible cross-correlations van- 
ishes when 7 — ► ±00 (See Ref . l25j) . Indeed, the term of 
Eq. equals to zero for e = y/2 — 1 « 0.42. 



0.2 ■ 



FIG. 4: Plot for the maximal Bell inequality violation B max as 
a function of the measurement time: f3 = t/tv- The analysis 
is performed in terms of current cross-correlations, the devia- 
tion from the coincidence of measurement times a = XB ~ XA is 



taken to be zero. The maximally possible violation is normal- 
ized to 1, thus the limit for Bell inequality violation is l/v2- 
The transmission probability e is the same as in Fig. 



ii) The case of short time measurements: f3 — * 0. 
Here we vary two parameters: e - the transparency of 
the superconducting fork, which defines the parameter 
7 = , — I — 5-, and a = XB ~ XA - which characterizes the 
lack of coincidence of the measurements by the two par- 
ties A and B. Again, the best violation occurs for a — > 0, 
and there is a dependence of maximal violation on trans- 
mission properties of the fork, which are described by 7. 
The Bell inequality in this case is the following: 



Bmax — 



2 a 



< 



1 

V5' 



(47) 



where we have chosen the Bell angles: Sip a b — Sip a b' = 
Sifia'b = 7r/4 — a, Sipa'v = 3tt / 4 — a, for which the viola- 
tion is maximal. The dependence of these angles on a is 
due to the presence of a under the cosine in Eq. (|46|l . 

In Fig. |21 we plot the dependence of the left hand side of 
the Bell inequality on a, for a certain 7 = 0.001. There 
is again a threshold in the transparency e ma x = 0.22, 
below which the violation is possible. For a transmission 
probability below the threshold e max the Bell inequality 

in Eq. (|47|l can be violated for: a < Ty\J (\/2 — l)/7 ~ 
0.6ry/e, while it remains inviolated in narrow regions 
separated by intervals: 

To compare the entanglement measure and the maxi- 
mal Bell inequality violation, we plot in the same Fig. 

It 



the concurrency C as well as the function y 1+ 2 C2 . 

is known [^ij that for mixed states B max is confined 
between these two values (for the case of pure states 

B m ax = \J 1+ 2°~ )> wn ich is clearly seen from the plots 
of Fig. According to our results, there may occur 
situations, where the state on the detectors is entan- 
gled (C ^= 0), but there is no Bell inequality violation 

{Bmax < 1/V2). 

From these two examples of i) a. = and ii) [3 = 
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FIG. 5: (Color online) Plot for the maximal Bell inequality 
violation B ma x depending both on the deviation from coin- 
cidence of current measurements: a = XB ~ XA and the mea- 
surement time /3r <C h/eV. The maximal value is normal- 
ized to unity. The transmission probability is taken to be: 
7 ~ e 2 = 0.001. Violation above the critical value of 1/V2 
occurs only at separate peak locations. 



we conclude that the violation of Bell inequality occurs 
for small measurement times: r < 2-ry and for mea- 
surements which are close to being coincident: \xa — 
Xb\/vf < 0.6ry/e. We performed numerical integrations 
of Eq. Ij46(l for arbitrary values of {a, (3, 7}. The results 
for 7 = 0.001 and 7 = 7 max ~ 0.1 arc shown on two- 
dimensional plots of Fig. [3] and of Fig. Note that the 
plot of Fig. represent a slice at a = of Fig. while 
Fig. 0] represents a slice at (3 = of the same figure. 
One notices that the oscillations of B max are damped 
as a is increased, and, as mentioned above that viola- 
tion is less likely when f3 is increased because the time 
window is too large. Such oscillations in the Bell pa- 
rameter were detected previously in the case of normal 
metal forks 01 ■ ^ is obvious from Fig. that, when the 
transparency of the superconducting fork is larger than 
the critical threshold, there is no violation of the Bell 
inequalities, although oscillations are still noticeable. 



VII. CONCLUSION. 

We have proposed a setup measuring energy entangle- 
ment for pairs of quasiparticles originating from single 
Cooper pairs in a superconductor. The use of energy fil- 
ters allows to separate the particles according to their en- 
ergy and to obtain orbitally entangled pairs of electrons, 
which arc then analyzed. This electronic setup and its 
detection apparatus is an analog of the momentum-phase 
entanglement performed in quantum optics |l5j . a) Two 
electrons from a Cooper pair correspond to the two pho- 
tons generated by down-conversion, which have different 
wave lengths. Indeed, Andreev reflection, or the emission 
of a Cooper pair, also involves two electrons with differ- 
ent energies, b) In optics the spatial separation of the 
photon beams follows automatically from the down con- 




FIG. 6: (Color online) Same as Fig. El except that the trans- 
mission probability is taken to be critical: 7 = y m ax ~ 0.1. 
The whole surface lies below the critical value 1/V2, so there 
is no violation of the Bell inequality. 



version process. Here, we used energy filters to perform 
such separation and considered short time measurements 
to allow interference of particles with slightly different 
energies. It is because the interference occurs between 
different energies that the "usual" zero-frequency analy- 
sis is not applicable here. This is why we developed a 
Bell inequality test for finite-times cross-correlations of 
electrical currents. 

First, we proposed a density matrix analysis for quasi- 
particles emitted from a superconductor and separated 
in energy. This allows to calculate the concurrence, one 
of the most widely used entanglement measures. Second, 
we showed how energy entanglement can be transformed 
into orbital entanglement and how it then can be tested 
via electrical current cross-correlations. For such cross- 
correlations, we constructed a Bell type inequality and 
showed what set of parameters lead to its violation. 

The noise cross-correlations between the left and the 
right reservoirs were derived in the context of the scatter- 
ing theory for normal metal-superconductor hybrid cir- 
cuits. Indeed, in stationary quantum transport these cor- 
relations enter explicitly the Bell inequality test. Here 
the manifestation of entanglement is found by adjust- 
ing the phases accumulated by electrons and holes in the 
conductor leading to the beam splitter (one can change 
magnetic fluxes within MZ interferometers by varying the 
magnetic field through them or by varying their area, 
like it was done in experiments with a single interferom- 
eter HU)- For a device consisting of perfect elements 
(filters, beam splitters,...), maximal Bell violation is ob- 
tained for short time of measurements: r < 2ry and 
small deviations from coincident measurements between 
the left and the right reservoirs: \xa~- xb\/vf < 0.6ry/e- 
for relatively small transparency e of the superconducting 
fork. We found that the violation of the Bell inequality 
depends on transmission properties of the fork, and it 
occurs for: e < e max « 0.22. 

A clear issue is to inquire whether such energy en- 
tanglement can be detected in experiments. Our setup 
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for generating entangled electron pairs has many compo- 
nents, each of which in principle can be built for instance 
using gated semiconductor hetero-structures. The design 
of a controllable normal metal-superconducting fork rep- 
resents a challenge but, recently, work in this direction 
has been promising yfj. All leads are assumed one di- 
mensional, and are assumed to have little or no backscat- 
tering as, for instance, in the experiments of Ref. |37j . 
The energy filtering fork can be constructed from an or- 
dinary fork with coherent quantum dots at two ends of 
it. The dots are adjusted to have transmission peaks 
symmetrically above and below the Fermi level. For bet- 
ter filtering these dots must have levels with an energy 
spacing larger than the bias voltage. 

Because of the complexity of this device, there are 
other complications which should be considered. The 
setup is composed of two MZ interferometers, which 
are both subject to dephasing jS^. The fluctuations in 
the gate voltages could trigger fluctuations in the path 
length, which consequently would add an uncertainty to 
the phase of electrons and holes, which enter the noise 
correlation expression of Eq. I|37|) . As long as the phase 
fluctuations are controlled (small compared to a phase 
angle 7r, meaning that the fluctuations in path length 



are below the Fermi wave length Xf), the Bell inequality 
will continue to be violated. Furthermore, a considerable 
precision is required in the building of such an interfer- 
ometer: the two arms need to have equal length, up to a 
distance vftv > Xp- 

In conclusion, we have proposed the first energy en- 
tanglement setup for testing the non-locality properties 
of electrons in nanostructures. With the previous work 
on spin entanglement using NS junctions, this latest work 
emphasizes the analogy with quantum optics. Our main 
goal here has been to show that energy entanglement can 
in principle be analyzed (concurrence) or detected (noise 
cross-correlation analysis) in this somewhat ideal device, 
using concepts and circuitry borrowed from mesoscopic 
physics. 
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APPENDIX A: DENSITY MATRIX 



Starting from Eqs. l(5jl- lfTT|) we average two anihilation operators products as: 

(c+^E^c^E^) - f m S mn S(E 1 - E 2 ) , 



(Al) 



where /,„ = fFD{E^eV) for electrons and holes, with /fd(-E) the Fermi-Dirac distribution. And averages of four 
operators are calculated using Wick's theorem: 



^ka[ { E x^i')c\ a , 2 {E 2 ,tji)c nia2 (E 3 ,t 3 )c niTl (E^U) 
= MEJfriE^S^S^JknS^SiEt - E 4 )S(E 2 - Eg) e -i*i(t.-t 4 0/» e +^ 8 (t,-V)/ft 
-fk(E 1 )f l {E 2 )5 <C72 8 (7 ^5 km 5i n 5{E l - E 3 )S(E 2 - £ 4 ) e -^fe-^)A e +^(V-**)A 

The resulting calculation leads to appearance of three terms with different spin symmetry: 

+ OG P + OO 



9(°>x) = 77 — / dExdE, 
(hv F ) 2 ^Ja Jo 



{<Wi < 



K, fc (x i 0^ fc (x i ))( % K^O^ i fe'))e-^ ltti - ti ' J e + ^^ t i'-*^/ fe (^i)/((^2) 



+(<, fe (* i 0^ fe (^)(^^ 

+(t**C Ci >S b ( a: 0)(«^^ 

+(«<' fc (x i 0^(xi))(< I (a; j )^ l (^))e +,Bl(t< - v) e- i ^ ( V-^)(l - /*(25i))(l - fi(E 2 )) 



(A2) 



(A3) 

(A4) 
(A5) 
(A6) 
(A7) 
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+(uJ fc (ziOuifc(*i))(«^ (A9) 

+(^ fe (^0^fe))(^(^O4K^'))e +iSl(tj - V) e +iE2( V^)(l - MBkMfr) (A10) 

+ (^ fc ( 2 : 4 0^fe))(^(^0^K^0)e + ^ l( ^-*-' ) e-^(V-*O ( ^ _ /fc (^))(l - / { (^))] (All) 

[Kk(*iO»i'*(ziO)M (A12) 

-K, fc (x^ fe feO)(^ (A13) 

-(^fc^O^fc^OJKK^KK^))^ 1 ^'"^ ^^-^^! - fk(Ex))fi(E 2 ) (A14) 

+(^ fc (a: < 0^fc(^0)W l (^)«i l (^))e^ l( V-V)e- iB »^- t ')(l - / fc (2Si))(l - /,(^))] } , (A15) 

where we have used the matix notation (cr^ = ±1): J„ 1(T2 = if y and fj, is a Pauli matrix. 

We consider the Andrecv reflection at the NS-boundary to be ideal (no normal reflection), and the reflection 
amplitude reads: ta = e - tarccos a. According to the properties of Bogolubov's equations the scattering amplitudes 
must satisfy the equations: 

S^(E) = Sf°*(-E), S%{E)=-S%*{-E), (A16) 
from this, one can derive the property: 

v ij/3 (x, E) = (-l) 5 ^ U y(x, -E), u ij0 (x, E) = (-l)^z^(x, -E), (A17) 
and using the identity fh(E) = 1 — f e (—E) one can simplify Eqs. i|A3(l - i|A15() : 

g{(J : X) ^(j^a'^a^a'^Gi'iG jjr S f7l a ^ S a2 ^ G i' j G iy ~\~ lo^o^l V a^a^Ei' j> E^ j : (A18) 



where 



1 

Gijixux^tj -U) = —Y dEf k (E)(u* k (x % )u 3k { Xj ))e- lE ^- u \ (Al9a) 

nv F ^ J-oo 
1 f +OC 

F ij (x i ,x j ,t j -t i ) = —Y^ / dEf k (E)(uU(x i )v jk (x j ))e- iE <- t >- t *l (A19b) 



k 



APPENDIX B: CURRENT CROSS-CORRELATIONS 



It is convenient to define the following matrix elements: 



A imin (E,E',x) = u in (E' \x)d x u* im {E,x) - d x u in (E' ,x)u* m (E,x), (Bla) 
x ) = v tn(E' , x)d x v im (E , x) - d x v* n (E', x)v im (E, t), (Bib) 
C imin (E,E',x) = u in (E',x)d x v im (E,x) - d x u in (E',x)v im (E,x). (Blc) 



These matrix elements have a useful symmetry: 



A imin {E,E',x) - —A* nim (E', E, x), (B2a) 
B imin {E,E',x) = —B* nim (E',E,x). (B2b) 



The answer for an average current is the following 

f + oo 



\ ie f +oa . 
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The calculation of irreducible cross-correlation gives: 
i- \ i- \ I » \ e 2 h 2 1 f + °° f +co 

(iiM)^^)) (U^U)) (ijfo,*;)) = 2 / dE dE'J2 (B4) 

' ^ ) Jo Jo m n 

f m {E){\ - f n (E'))e^ E - E 'X t '- t *V H [A im i n (E,E',x i ) + Bt min (E,E', Xi )] [A* mjn (E,E', Xj ) + B jmjn (E,E', Xj )] 

+f m (E)f n (Eyr« E+E '^-^/ h C; mm (E',E,x t ) [C jnjm (E',E, Xj ) + C jmjn (E,E', Xj )} 

+ (1 - f m (E))(l - UE'))e + ^ E+El ^-^/ h C imm {E,E\ Xl ) [C* mjn (E, E' , Xj ) + C* njm (E' , E, Xj )}~ . 

The averages of products for two and four anihilation operators were calculated using Wick's theorem in Eq. I)A2J) . 
In Eq. i|B4|) all terms like in Eqs. (|Bla|) - (|Bla|l may be written in terms of the scattering matrix elements 

Aika,up{E, E , x) 

Bika,iip{E, E , x) 
Cika,il@(E, E' , x) 

here we have linearized the k- vectors in energy, thus we are not interested in the quadratic behavior of the spectrum 
near Fermi level and hence we neglect the dispersion of elec tron and hole wave packets. 



= 2ik F st%*{E)sf(E')e +i{B - B ' ) ^ - 2ik F 5 ae 5p e 5 lk 5 a e- l{E - E,) ^ (B5a) 
= 2ik FS h l £{E)sf*{E , )e- l(E - E ' ) ^ - 2ik F 6 ah 6p h 6 ik 6 il e +i{E - E ' ) ^ (B5b) 
= 2ik F s^(E)sf {E')e- i{E+E,) ^ - 2ik F 8 ah 5 Pe 8 lk 8 a e +l(E+E ' ) ^ , (B5c) 
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